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Abstract 

A quantum mechanical treatment of an asymmetric double-well potential 
(DWP) interacting with a heat bath is presented for circumstances where the 
contribution of higher vibrational levels to the relaxation dynamics cannot be 
excluded from consideration. The deep quantum limit characterized by a dis- 
crete energy spectrum near the barrier top is considered. The investigation is 
motivated by simulations on a computer glass which show that the considered 
parameter regime is "typical" for DWPs being responsible for the relaxation 
peak of sound absorption in glasses. Relaxation dynamics resembling the 
spatial- and energy-diffusion-controlled limit of the classical Kramers' prob- 
lem, and Arrhenius-like behavior is found under specific conditions. 
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I. INTRODUCTION 



Relaxational dynamics of materials which evolve under slow structural change at low 
temperature can frequently be studied in a double-well potential (DWP) picture that allows 
tunneling between the two wells. A prototype of such a system, important in many chemical 
and biological systems, is the transfer of hydrogen atoms along preexisting hydogen bonds 
for instance in crystalline benzoic acid and carboxyl dimers.B~i Another example is the 
tunneling of Si0 4 - or Ge0 4 -tetraeder in amorphous Si0 2 and GeO 2 .0i 

One of the key issues in understanding the relaxation process is the coupling of the motion 
of the light tunneling particle to the havier surrounding atoms constituting the heat bath. 
In a condensed phase the coupling leads to structural rearrangements of the environment 
which provides a mechanism of relaxation for the system. 

The theoretical description of dissipative DWP-dynamics has been widely elaborated in 
the low temperature limit where a two- level description (spin-boson model) is sufficient. 
At higher temperatures where the excitation of vibrational levels can no longer be neglected, 
the effect of intra-well relaxation has been given explicitly by Silbey and coworkers^ for 
a symmetric DWP. Lateron Meyer and ErnstH considered a biased DWP with crystalline 
benzoic acid dimer as an example. In both papers the investigations have been restricted 
to a regime where the tunneling splitting A n of all relevant doublets is still very much less 
than the interdoublet spacing, i.e., A n <C uj where ujq is the frequency of small oscillations 
around the two minima of the DWP. They found in the limit of fast intra-well vibrational 
relaxation - where the population of the doublets is in thermal equilibrium - that the escape 
rate is the thermally averaged tunneling rate into the lower well, 



In many practical applications, however, whenever the two-state approximation seemed no 
longer sufficient, a simple classical relaxation over the barrier V with an Arrhenius rate 
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has been used successfully.!!!'! 

Hence, concerning the understanding of relaxation dynamics in many chemical and phys- 
ical systems at higher temperatures, two questions are of major importance: (i) How dense 
are the energy levels of the DWP as compared to the potential height, i.e., is the classi- 
cal picture underlying Eq. ( |1.2[ ) justified, or does the relaxation dynamics depend on the 
quantum mechanical eigenvalues of the DWP? (ii) Is the DWP coupled to bath modes with 
frequencies of the order of its potential height, and does the relaxation dynamics depend on 
the parameters of the heat bath? 

In this work these questions are analyzed for the specific example of DWPs in a structural 
glass. There DWPs correspond to the motion of a local group of atoms or molecules between 
two local energy minima. Low barrier DWPs are known to be the source for the tunneling 
properties observed in glasses like SiC>2 below 10 K.0 In contrast, DWPs with higher barriers 
(in SiC^: (V^, e ak) /&b ~ 500 K) are believed to be responsible for the relaxation peak of sound 
absorption, observed around 50 K in SiC^.Bi Our goal is to answer the above questions for 
these DWPs in order to obtain a closer understanding of the relaxation dynamics at the 
absorption peak. 

To this aim we first perform in Sec. 2 computer simulations on a model glass to de- 
termine the eigenvalue spectrum and the interaction with phonons of those DWPs which 
are responsible for the relaxation dynamics in glasses at higher temperatures. We find that 
DWPs relevant at the sound absorption peak only contain a few energy levels below the their 
barrier height V and are coupled to bath modes with frequencies as large as V. This serves 
us as a motivation to analyse the relaxation dynamics in a single DWP as depicted in Fig. 
1 which is characterized by the parameter regime 



k B T ~ Uujq ~ V , (1.3) 

for instance, T « 30 — 50 K, V/k-Q « 300 K and V/Uujq « 2. In Sec. 3 a simple model 
is presented which allows an analytical solution of the relaxation dynamics by applying 
standard approximation, as shown in Sec. 4. With this model we can discuss in Sec. 5 the 



apparent paradox that a classical Arrhenius rate does well in many practical applications 
although the relaxation dynamics is far from being classical. In Sec. 6 we close with a short 
summary. 



II. SIMULATIONS 

In previous work Heuer and Silbey have developed an algorithm which systematically 
searches for DWPs in glasses simulated on the computer .ElIE This algorithm has been readily 
applied to the specific case of NiP which can be described as a binary Lennard- Jones-type 
glass. For general reasons one expects that Vp ak is approximately proportional to the glass 
transition temperature.@ Hence from T 5 (NiP)/T 5 (Si02) ~ 2/3 one may roughly estimate 
Vpeak (NiP)/&B ~ 300 K. We proceed in two steps. First we present a rough estimation of 
the number of energy levels in DWPs with V ~ V^ ea k- Second we explicitly determine for 
the NiP computer glass the interaction of the DWP with the heat bath. 

For the estimation of the number of energy levels below the potential height it is sufficient 
to consider a symmetric DWP. The potential energy along the reaction coordinate x, which is 
defined as the collective coordinate connecting the two local energy minima, can be expressed 
quartic polynomial 

V{x)= l -^{x-d/2)\x + d/2f (2.1) 

where V is the potential height and d the distance of both minima. We chose the DWP 
symmetric around x = 0. The length scale d is defined such that the collective dynamics 
along this DWP can be viewed as the dynamics of a single particle with mass m where m is 
the average mass of the atoms (for NiP m = 56m p with m p the proton mass). 
It is straightforward to calculate the vibrational frequency ujq in both wells as 

= < 2 - 2 ' 

In Ref .0 it has been shown that for NiP the value of the prefactor 16V/ d A of the quartic term 
for symmetric DWPs on average is given by BA^/a? with B/ks = 4.3 x 10 4 K, A4 = 0.39 



and a = 2.2 x 10 10 m. Elimination of d in Eq. (|2~^ ) yields uj^ = 8^VBA 4 /(ma 2 ). Inserting 
the numbers for NiP given above yields for V/k-Q = 300 K that 

V/hu « 5.5 . (2.3) 

This value should be viewed as an upper limit since for DWPs with large barrier heights the 
quartic contribution is expected to be even larger than the average contribution. However, 
this would only reduce the number of energy levels below the barrier height. 

For evaluating the coupling to the bath modes we analyzed a NiP computer glass with 
iV = 500 atoms using periodic boundary conditions. The simulations have been performed 
at zero temperature; details can be found in Refs. |T^JT6l . To first approximation the energy 



surface of the computer glass adjacent to the DWP can be written as 

H = H (x) + xJ2 \V 3 + \m E u)y) (2.4) 

3 3 

where H (x) describes the relaxation mode in the DWP, Uj the frequencies, jjj the coordi- 
nates of the bath modes, and Xj the coupling of the bath modes to the DWP. This potential 
energy term contains the full anharmonicity with respect to the DWP mode. The other de- 
grees of freedom are only considered up to their harmonic contributions. Of course, in order 
to obtain the full Hamiltonian of the system one has to add the kinetic energy of the relax- 
ation mode and of the bath modes. For our simulation we chose a DWP with V/k-Q ~ 300 K. 
Starting from one minimum we calculated all second derivatives of the energy with respect 
to the positions of all 500 atoms, yielding a 3iV x 3iV dimensional matrix A. Let d denote 
the unit vector in the 3iV-dimensional configuration space connecting both minima of the 
DWP. The matrix A is diagonalized in the space orthogonal to d , yielding eigenvalues cij and 



eigenvectors ej. The corresponding eigenfrequencies Uj are given by luj = ^aj/m. Three 
eigenvalues equal to zero reflect the periodic boundary conditions. The coupling parameters 
Xj are calculated via Xj = dAij. It is worth noting that our vector d which is just a straight 
line path between the two minima of the DWP is not the reaction path between the two 
wells, which would be curvilinear. However, the simulations in Refs. [H)| , |TE| have shown that 
the curvature of the reaction path is not very large. 



Later on it turns out that the relevant quantities describing the coupling strength of the 
bath modes to the DWP is the spectral densities 



J(ui) is plotted in Fig. 2. One can clearly see that bath modes exist for all frequencies 
smaller than 400 K. As a result, DWPs relevant at the sound absorption peak which only 
contain a few energy levels below the their barrier height V are coupled to bath modes with 
frequencies as large as V. The high-frequency modes are somewhat more localized that the 
average modes and contain the dynamics of the order of 100 atoms. Such an effect has been 
already observed in previous simulations.0 

A closer inspection of the data shows that the coupling is strongly anharmonic. This 
can be easily checked by repeating the above procedure calculating the second derivatives 
at different points along the reaction coordinate. It turns out that the A? and Uj of the 
individual modes strongly vary along d whereas the distribution function J(u) is insensitive. 
The reason for this is that J(oo) arises from an ensemble average over a large number of 
oscillators whose parameter variations are statistically independent in a first approximation. 
However, one should keep in mind that due to the finite size of the simulation box we 
do not obtain the low-frequency phonons in the Kelvin-regime which are relevant for the 
transitions between both minima at low temperatures. For those phonons one expects that 
the anharmonic contributions are significantly smaller. A more detailed discussion of the 
anharmonicities is beyond the scope of the present paper. 



(2.5) 
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where 




(2.6) 



III. THE MODEL 



We now construct a simple model which governs the relaxation dynamics of the relevant 
DWPs. For reasons of simplicity we restrict ourselves to only two pairs of tunneling doublets 
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below the barrier. 

The phonons couple to the tunneling coordinate. Their effect is two-fold: They induce 
vibrational transitions between the eigenstates of the DWP and destruct the tunneling 
coherence within each doublet. Concerning the first effect we assume that barrier penetrating 
transitions are negligible. The relevant matrix element between the localized states reads 



M = (Oa\H\la) = Y, c f( b i + h )) t 3 - 1 ) 
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with phonon operators satisfying = 6jj> and the coupling constants Cj = Cj of each 

mode. Here, the first index labels the doublet and a =L, R the left and right well. The 
second effect of the phonon coupling is modeled by a diagonal coupling which detunes \nL) 
against |nR) (n,m = 0, 1). The corresponding matrix elements read 

e = l({nR\H\nR) - (nh\H\nh)) 

= £4% + &]) (3.2) 

j 

with (note that |((nR|x|nR) — (nL\x\nL)) « d) 

= k . (3.3) 



(0a|x|la) 



C W " d 



Corresponding to Cj and cp we also define the spectral densities J v (u>) = J(u>) and J t (u) = 
^J(w) where J(u) has been defined in (|2.5|) . 

In case of our computer glass the factor k can be estimated for symmetric DWPs as given 



by Eq. (|2.1| ). Since (0a|x|la) = Jh/2rruvo one directly obtains from Eq. (|2.2| ) and Eq. ( p.3|) 
k ~ 0.05. Hence, though the intra-well vibrational transition M may be considered as a 
weak perturbation, the diagonal coupling e, generally, may not. The main reason for this is 
the "polaron effect", i.e., the slowing down of the tunneling transition associated with the 
shift of the oscillator coordinates into the other well. 

As mentioned in Refs. a suitable decomposition of the Hamiltonian for second order 
perturbation theory is achieved by the polaron transformation!!! (er™ = |nL)(nL| — |nR)(nR|) 
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(3.4) 



We denote the energy separation of the pair of tunneling doublets bylll AE = huo m /2 
and the asymmetry by A < hoooi (otherwise the states isolated in the lower well will drop 
out of the dynamics and the doublets have to be redefined). Now we may eventually define 
our model Hamiltonian. With the shift operator 



2c (t) 

B ± = exp ^±E^-( fe i- fe 5) 



(3.5) 



and the tunneling frequency renormalized by a Debye- Waller factor (-B±)b = e 
Appendix A), 
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where 5B± = B± — (B±) B , and H B = Y^j^jbjbj. Our goal is to apply second order 
perturbation theory in H mt . In this approximation the intra- well dynamics is governed by 
the one-phonon transition between the vibrational levels, i.e., by the Orbach process. In case 
of more than two doublets below V, we have to sum over all pairs (see below). We note in 
passing that the spectral density J(u>) only holds for transitions from level to level 1. For 



transitions changing the vibrational quantum number more than one the transition matrix 
elements are only due to the anharmonicity of the potential and the interaction. This may 
lead to a somewhat decreased spectral density. However, since we are only interested in an 
order of magnitude estimation of the coupling to the bath we do not analyse the influence of 
anharmonicies on the coupling in greater detail. Rather we state that for the observed weak 
intra-well coupling (as compared to inter-well coupling due to k <C 1) the Orbach process is 
indeed the dominant relaxation mechanism. 



IV. DYNAMICS 

To determine the time evolution of the tunneling doublets we make use of Mori's pro- 
jection operator technique^. We adopt throughout the paper the fast vibrational relaxation 
limit which is relevant for most experimental situations. Denoting the downwards vibra- 
tional transition rate between the two doublets by 701 and the tunneling escape rate in the 
other well (inter-doublet transition rate) by i~o and I\, this means that we assume 

A, A < 701 • (4.1) 

This guarantees local thermal equilibrium in each well. The relaxation rates are defined 
in Appendix A. For weak coupling between the vibrational levels and in the temperature 
regime ( |1.3|) , the upwards vibrational rate 710 = 701 e' 13 ^ 01 is typically much smaller than 
the level spacing 

e n = ^Ayh 2 + Al (4.2) 

of the two doublets [cf. Fig. 1]. Furthermore, due to the strong increase of the tunneling 
matrix element near the barrier top, the ground state spacing is much less than the spacing 
near the barrier top. Hence we are interested in the regime 710 <C eo <C e\ with no fixed 
relation between 7 01 and e , e\- 

Relaxation dynamics is probed in long time or low frequency experiments where the time 
scale of the experiment is far off any resonant time scale of the system. Thus, any coherency 
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has already decayed and the system shows pure decay. The experimentally accessible quan- 
tity is the population difference between the two wells. The relevant operator measuring this 
quantity is 

(I \ 
0-100 
10 
\0 0-1/ 



Q 



(4.3) 



Its thermal expectation value reads 



(Q) = -tanh(/M/2) . 



(4.4) 



In linear response theory all information is contained in the symmetrized correlation function 
of the relevant operator 



C(t) = \(5Q(t) 5Q) + t^-t 
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gSQe~ iCt 5Q 
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(4.5) 



where 5Q = Q — (Q), and £ = U~ [H, *] is the Liouvillian. Due to local thermal equilibrium 
in each well, the density matrix is given in the canonical form g = exp(—j3Hd) /Tr exp(— (3Hd) 
with a Hamiltonian Hd which contains only the diagonal part of Hq in the localized basis. 
The spectral function C(u>) = (1/2) /f^ C(t)e luJt dt can be directly measured in neutron scat- 
tering experiments, or via the dynamical susceptibility in acoustic or dielectric experiments. 

In Appendix B we calculate in second order perturbation theory in H [nt the correlation 
function (|4.5f) by using Mori's continued fraction representation for the complex correlation 
function C(A) = Jq 00 e~ xt C(t)dt (Re(A) > 0). After inverse Laplace transformation we find 
in the limit A — > [cf. Eq. B.5] our final result 



C(t) = ((SQ) 2 ) e 

with ((SQ) 2 ) = sech 2 (/3Z\/2) and the escape rate 

k — kt -\- k v . 
9 



-2kt 



(4.6) 



(4.7) 



The tunneling rate k t is given by the usual small polaron expression averaged over the 
thermal level occupation 

h = i(r + r x ) p 

^ l°° dt (A 2 + A*)? cos(At/h) (cos^(t)] e»« - l) (4.8) 



where 



< 2 «>" = tS^ (4 ' 9) 



In the fast vibrational relaxation ( |4.1j ) the effect of the vibrational transitions enters k t only 
through the thermal occupation number of the localized states. 
The vibrational relaxation rate 



j4f 7oi_ 

2 4 + 7oi ^ hW(n + 1 



is characterized by the ratio of the doublet splitting ei and the equilibration rate 701 between 
vibrational levels. This type of expression for k v has first been proposed by Sussmannil on 
ground of a simple wave function argument. For A — it is identical to the transition rate 
between E and A symmetric states in a rotational tunneling problem.Bi! 



V. DISCUSSION 

A. Analogy to the spatial- and energy-diffusion-controlled limit 

Our results can be summarized in the following physical picture. There exist two relax- 
ation mechanism for a particle in the energetically instable well : (1) Incoherent tunneling 
with rate kt within the thermally occupied tunneling doublets, or (2) vibrational transition 
to a non-thermally occupied doublet of higher energy, coherent tunneling within the doublet 
and vibrational decay to the lower well which results in a rate k v . This interpretation also 
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follows nicely from the projection operator method used in Appendix B. Relevant degrees 
of freedom are the diagonal terms in H , i.e., the population in the left and right well. In a 
projection operator formalism, relevant degrees of freedom (projected out by V at an earlier 
time t' < t) couple via an interaction (£) to irrelevant degrees of freedom (projected out by 
<2), which subsequently evolve in time (Q£Q) and, due to a second interaction (£), acquire 
relevancy again, thus influencing the evolution of the relevant degrees of freedom at the 
present time t. Eqs. (B.6) and (B.7) are just the Laplace transform of the kernel describ- 
ing this process. Now, k t propagates that part of the irrelevant degrees of freedom which 
destruct coherence of the relevant degrees of freedom, i.e., A n 5B±, whereas k v propagates 
that part which maintains coherence, i.e., A n . 

Process (1) has been discussed by Parris and SilbeyEH and Meyer and Ernst!. This 
process always prevails if 701 3> e% because then the particle has no time for tunneling 
before decaying from the upper doublet. In this case the preexponential factor of k v 



controlled limit in the classical Kramers' problem.E3 The escape to the other well is hindered 
by the damping of the tunneling motion, i.e., k ~ k t . The tunneling escape rate k t is given 
by the Boltzmann-weighted average of the tunneling equilibration rate within each tunneling 
doublet at given vibrational level. This equilibration is hindered by the distortion of the 
lattice in the tunneling process, i.e., by the polaron effect. The suppression of the vibrational 
rate for 701 ^> ei can be thought off as a phonon bottleneck effectE The physical reason for 
the suppression of k v is that phonons resonating between the upper and the lower level of 
two different tunneling doublets are not distinct because the level width exceeds the doublet 
spacing. As a result one cannot add up their rates. Instead they can interfere resulting in a 
suppression of k v . 

Process (2) has first been proposed by SussmannS on ground of a simple wave function 
argument. It prevails provided that (z) the vibrational transition to the upper doublet is 



foi — 




(5.1) 



is approximately / 01 



A* 
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faster than the direct decay via incoherent tunneling, 



7io = 7oi e 'bt > r c 



(5.2) 



and, necessarily, if (ii) the decay of the excited state is slower than the tunneling process, 
7oi < E\. In this case the particle has sufficient time for oscillating coherently fro and back 
in the upper level, and, therefore, a maximal, i.e., 50% probability for decaying in the other 
well. As a result the energy transport to the barrier top becomes the rate limiting process, 
i.e., k ~ k v . This is the situation analogous to the energy- diffusion- controlled limit in the 
classical Kramers' problem. It should be noted that the quantity which determines the 



frequency of coherent oscillations is the doublet spacing e\ = \JA\ + A 2 /Ti 2 and not the 
tunneling frequency A\. Hence, the asymmetry acts in suppressing the bottleneck effect. 

Clearly process (2) becomes increasingly likely with increasing temperature. Eq. ( |5.2| ) 
defines a transition temperature T*. Furthermore a large e\ is needed. Hence, upper doublets 
with energies of 0(V) are preferable. 



Now we present a derivation of an apparent Arrhenius behavior at higher temperatures 
with an activation energy 0(V). We still consider the situation where the eigenvalues are 
discrete and not dense at E m V. Our previous results are now generalized to many pairs by 
noting that the total rate is the sum over all pair rates weighted by the Boltzmann occupation 
factor of the initial state. We denote the (downwards) equilibration rate between the mth 
and the nth vibrational state (n < m) by 7 nm . This provides the rate 




B. Apparent Arrhenius behavior 



k = k t + k, 
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with 



nm 




(5.4) 
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For [n = 0, m = 1) the term for k v boils down to Eq. ( [4.11] ) . Now we consider the relaxation 
of a single DWP with given parameters. As discussed above for T > T* the relaxation 
is mainly determined by k v . In order to obtain the temperature dependence one has to 
check by which doublets (n, m) this term is dominated. As a first approximation we neglect 
the dependence of 7„ m on n and m, hence introducing a single value of 7. Then f nm only 
depends on m. This may be a poor approximation. However, since other quantities like the 
tunneling matrix element or the Boltzmann factor exponentially depend on the indices this 
approximation is sufficient for our present purposes. 

We now ask by which value of m the term f nm is dominated. We first note that only 
doublets with e m > 7 can avoid the bottleneck effect, and that, even among those, doublets 
with A m j A <C 1 are unimportant. Evidently, f nm cannot grow beyond 7/2 — the strict 
maximum is only reached for A 2 m — > 00. However, already for 

^ = ^ 2 A 2 + 7 2 (5-5) 



the value 7/4 is reached. Hence, Eq. (|5.5| ) is an appropriate criterion that f nm is close to 
its maximum value. One can estimate from the numerical data that the average value of 
it J (uj)k-B / Ti is of the order of 30 K. Having in mind that this estimation may be somewhat 
too large (see Sec. 3) we may estimate ^h/k-Q = O(10 K). This is somewhat smaller but of 
the same order as the numerical value we have found for Uq. Hence, at least for asymmetric 
DWPs, criterion ( |5.5|) can hardly be fulfilled for real systems. However, doublets which come 



closest to (5.5) will still dominate f nm . Evidently, these are levels m close to the barrier. 
This is related to the exponential dependence of the tunneling matrix element on the energy 
of the corresponding level. Only close to the barrier the maximum value A ~ 2ujq/tt is 
reached.0 

Additionally taking into account the Boltzmann factor one can see that all levels n 
contribute as long as exp(/3E m ) 3> exj>(/3E n ). For (3E m » 1 this implies that most levels 
n indeed contribute. 

Summarizing the arguments presented above we obtain 
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k v « r, 1 e~^ u (5.6) 

with an activation energy AU = 0(V) and a prefactor Tq 1 which is of the order 7/4 x m. 
The latter factor results from the summation of most of the states below the m-th level. 
Since our numerical results indicate that 7 is somewhat smaller but close to Uq we may write 
Tq 1 m 0(u!o/2tt). As a surprising result, for sufficiently high temperatures, T > T* such 
that k v 3> kt, the relaxation rate between the two wells is an effective Arrhenius rate where 
both the prefactor and the activation energy are of the same order of magnitude as in the 
classical limit. 



VI. SUMMARY AND CONCLUSION 

The intension of this work was to investigate the remarkable success of fitting relaxation 
data of disordered solids at higher temperatures by a simple Arrhenius rate despite the fact 
the underlying classical picture might not be justified automatically. 

Due to the lack of knowledge of microscopic DWPs in disordered solids we have used 
numerical simulations on a computer glass to characterize the DWP and the bath modes 
before solving the dynamics for this particular situation. It should be noted that these sim- 
ulations have proven to be reliable by reproducing characteristic low-temperature properties 
of glasses. 00 Hence, the parameter regime for which we have solved the dynamics can be 
considered as "typical" for structurally disordered solids. The simulations have taught us 
that DWPs relevant at the sound absorption peak, i.e., those with large V, only contain a few 
energy levels below their barrier height V and are coupled to bath modes with frequencies 
as large as V. 

To resolve the paradox to the apparent classical behavior we have constructed a simple 
model of two tunneling doublets which are coupled by intra-well vibrational transitions. 
Independently of how closely this model reflects the true physical situation in detail, it has 
provided a simple physical picture how even in this deep quantum mechanical regime an 
apparent Arrhenius behavior might emerge. Solving our model by standard approximations 
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we could confirm the suggestion of SussmannE± at high temperatures the vibrational transi- 
tion to a doublet near the barrier top, followed by coherent tunneling between the wells and 
vibrational decay to the bottom becomes a faster relaxation mechanism than direct decay 
via incoherent tunneling. Clearly the derivation which eventually led to the Arrhenius rate 
( |5.6p is very scetchy. In general there will be corrections to the term satisfying ( |5.5| ) which 



are in the present formulation difficult to quantify. Furthermore, if more than two doublets 
are involved, it is difficult to define quantitatively the transition temperature T* which sep- 
arates the tunneling from the vibrational dominated regime. Hence, generally, both terms 
k t and k v may contribute and the Arrhenius rate has to be viewed as an idealization. 

However, we believe that the picture drawn provides a physical intuition on the remark- 
able fact that a simple Arrhenius rate is so successfull in explaining relaxation data even 
just above temperatures where the quantum mechanical tunneling effect had dominated the 
relaxation dynamics. In addition to this it is reasonable to expected that in disordered solids 
like glasses the detailed form of the barrier distribution function has a much greater impact 
than possible corrections to the Arrhenius rate. Hence, the customary application of this 
rate for practical purposes can be justified a posteriori by the present picture even in the 
deep quantum regime. 
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APPENDIX A: BATH CORRELATION FUNCTIONS AND RELAXATION 

RATES 



Assuming that the bath is not disturbed by the DWP we can calculate bath correlation 
functions with Hb and qb = exp(— (3Hb)/Zb- Defining 
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g{t) = f dcu J^w) (n(cu)e^ + (n(w) + l^"*"*) (A.l) 

^ (f) = 2 f° dcu^M sin(wt) (A.2) 
Jo u 2 

roo J-t(uj) 

g 2 (t) = 2 du;-\±coth(f3hw/2)cos(u;t) (A.3) 
Jo UJ 2 

W(T) = y 2 (0) (A.4) 



the bath correlation functions read 



{B + (t)B_) B = (B_(t)B + ) B = (B + B_(t)) B = (B_B + (t))\ 



= exp 



g 2 (t)-i gi (t)}e- 2W ^ (A.5) 



(M(t)M) B = g{t) (A.6) 

Furthermore, we need the following modification of the vibrational relaxation function 

g ± (u Q , t) = e^ ' f J™ du J v (u) (n{u)e ±iu)t + (n(u) + l)e^*) . (A.7) 

In the following two kind of relaxation rates are needed. First, the transition rate within the 
nth doublet (inter-well transition) 

r n = ^ J°^dt ~A\ co S (At/h) (cosM0]e 92 W - l) . (A.8) 

This expression is well-known from small polaron theory. 

The second rate is the transition rate between the doublets (intra-well transition) 

7io= 2^/" dt(g_(u 01 ,t) + g + (-u 01 ,tj) (A.9) 
7oi = 2F J™ d*( ff _(-w i, *) + g+(woi, t)) . (A.10) 

7io is the transition rate from the zeroth to the first vibrational level and given by 

7io = 7rJ„(a;oi)n(a;oi) . (A- 11 ) 

Here, n(uoi) = [exp( / 9^a;oi) — is the Bose factor. 701 and 710 satisfy the principle of 
detailed balance 710 = 701 e _/3n ^ 01 . 
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APPENDIX B: FORMULATION OF THE DYNAMICS IN THE MORI SCHEME 



If we define the symmetrized scalar product (A\B) := \{A^B + BA^) we can write 
the complex correlation function C(A) = / °° e~ A 'C(t)dt (Re(A) > 0) as a resolvent matrix 
element 

C(X) = (5Q\[X + tC]- 1 \5Q) . (B.l) 

We define the projector V = \5Q)ri~ 1 (5Q\ = X — Q and its complement Q with the normal- 
ization 7] = (SQ\SQ) = sech 2 (/3Z\/2). The resolvent identity 

[A + iC VP + C PQ [X + iCqq^Cqv] V[X + iL\- x V = V (B.2) 

where C AB = ACS with A,B G {V, Q} defines a relaxation kernel fc(A). With (5Q|£|5Q) = 
we find 

C < A > - j^hxj < a3 > 

where the relaxation kernel is given by 

2k(X) = (5Q\C PQ [X + iC QQ ]- l \C Qr 5Q) v' 1 • (B.4) 

For long times the transition rate k between the left and right well is given by the limit 
A -> 

k = limA;(A) . (B.5) 

To calculate this function, we first separate k(X) into two terms k(X) = k t (X) + k v (X) with 

2k t (X) = (5Q\C VQ [X + iCqqY^R) rf 1 (B.6) 
2k v (X) = (SQljCrQiX + ijCQQ^lRjr]- 1 (B.7) 

where 5R contains all matrix elements of Cq-pSQ arising from H® t and R contains all matrix 
elements of Cqj>8Q arising from the non-diagonal part of H . Note that does not 

contribute to Cg-pSQ. 

17 



We now apply lowest order perturbation theory in H int = + H^, i.e., in A n SB± 
and M. Since terms which are proportional to the product of these two perturbations are 
of higher order, we neglect the vibrational coupling M in the tunneling part k t , and the 
tunneling coupling A n 5B± in the vibrational part k v . 

It is a standard result of spin-boson literature that for a single tunneling doublet the 
tunneling transition rate defined by k t = k t (X = 0) is given by expression (A.8) in second 



order perturbation theory in the tunneling matrix element. Repeating the same kind of 
calculation for the pair of doublets one easily finds (|4.8|) in case of fast vibrational relaxation 
(0). 



The calculation of the vibrational rate is essentially identical to the calculation of the 
transition rate between E and A symmetric states in a XH 2 rotational tunneling problem. 
Let us introduce the notation that E nm denotes a matrix with a one in the n-th row and the 
ra-th column and zeros elsewhere. With this, we separate R into two part R = Zio-Ro + ^i-Ri 
with R = E 2 i — E%2 and R\ = E i3 — E u . Inserting this and the projector 

V'= J2 \Ra)Va\Ra\ (B.8) 
a=0,l 

with T] = [1 + e"^ 01 ]" 1 and 771 = [e^ 01 + l]^ 1 into k v (X) and noting that 
(5Q\£\R )r]o 1 ri~ 1 = A and (SQ\C\Ri)rji 1 ?7 _1 = A 1 one finds that the vibrational part 
of the relaxation kernel becomes 

fc«(A)«2?Gii(A) (B.9) 

with 

G n (X) = (R^X + iCqq}- 1 ^) . (B.10) 

Because of A\ ^> A , we have neglect all terms oc A . The calculation of this function 
in second order perturbation theory in the vibrational coupling M allows the replacement 
Cqq — > C = H^IHq + Hfol, *] where H + H^t have been defined in (|3.7| ). The problem 
is now essentially idential to the calculation of the scattering function in a XH2 rotational 



tunneling system. Details of this can be found in Refs. |24|j25| . Following essentially their 
lines one finds with the definition k v = k v (X — 0) the result ( |4.11| ). 
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FIGURE CAPTIONS 

FIG. 1: Schematic representation of a biased double-well potential with V ~ 2hu . The ground 
state and first excited state tunneling doublets are separated by an energy AE hu . 

FIG. 2: Spectral density J v (oo) = J(oo) for the NiP Lennard- Jones computer glass. 
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